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, Abstract 



To perform a queuing analysis or design in a communications context, we need to estimate 
the values of the input parameters, specifically the mean of the arrival rate and service time. In 

■ this paper, we propose an approach for estimating the arrival rate of Poisson processes and the 
average service time for servers under the assumption that the service time is exponential. In 
particular, we derive sample size (i.e., the number of i.i.d. observations) required to obtain an 
estimate satisfying a pre-specified relative accuracy with a given confidence level. A remarkable 
feature of this approach is that no a priori information about the parameter is needed. In 

■ contrast to conventional methods such as, standard error estimation and confidence interval 
0^ , construction, which only provides post-experimental evaluations of the estimate, this approach 

\ allows experimenters to rigorously control the error of estimation. 
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The assumption of Poisson arrivals is usually valid in modeling traffic in communications pQ. 
There is well-developed mathematical theory justifying such an assumption. Also, it is often 
assumed that the service time of a server in a queuing system is exponential. In the case of 
telephone traffic, the service time is the time for which a subscriber engages the equipment of 
interest. In a packet-switching network, the service time is the transmission time and is therefore 
proportional to the packet length. It is difficult to give a sound theoretical reason why service 
times should be exponential, but the fact in most cases they are very nearly exponential [T]. For 
design and analysis purpose, it is desirable to estimate accurately the load currently generated 
by each device inter-connected in a network. For example, the rate of packets generated by a 
terminal or the size of packets. 



*The author is currently with Department of Electrical Engineering, Louisiana State University at Baton Rouge, 
LA 70803, USA, and Department of Electrical Engineering, Southern University and A&M College, Baton Rouge, 
LA 70813, USA; Email: chenxinjia@gmail.com 



1 



2 Error Control 



Conventional advice recommends approximating the Poisson arrival rate A as S where n is the 
number of arrivals during an interval of time with length T. The problem with this method is 
that we do not know how accurate the estimate is. Moreover, we have difficulty in choosing the 
length T to guarantee a certain accuracy. Observing over-long time interval is wasteful. On the 
other side, inadequate observation will lead to poor estimate. 

Situations are similar for estimating the average service time \i. The estimate of (i is calculated 



where Xi, i = 1, ■•• , n are i.i.d. observations of service time. For the purpose of accuracy 
evaluation, the standard deviation of /2 is estimated as 



It should be noted that such an estimate of standard deviation is not a good measure of accuracy, 
since itself is also a random variable. 

An alternative of accuracy evaluation is confidence interval construction. The limitation of 
standard deviation estimation and confidence interval construction is that they are only post- 
experimental evaluation of accuracy based on fixed number of observations. These two conven- 
tional methods do not control the accuracy of estimation because the error is determined once 
the number of experiments is fixed. There are two natural methods of error control. First, one 
can specify an absolute error bound e > and confidence parameter 5 S (0, 1) and determine the 
sample size n such that 



where rj is an estimate of the parameter rj (e.g., Poisson arrival rate or average service time). 
Second, one can specify a relative error bound e G (0, 1) and confidence parameter 5 G (0, 1) and 
determine the sample size n such that 



The problem with the first method is that it is not feasible to determine the sample size without 
a priori information of the parameter. Moreover, an absolute error bound is not a good indicator 
of the precision of the estimate. The second method is better since a perfect measure would be 
expressed in terms of a relative error bound. More importantly, the sample size can be determined 
without any information of the parameter to be estimated. 

3 Arrival Rate 

To estimate the arrival rate of a Poisson process, we have 




Pr \\r\ - 7/| < e} > 1 - 5 
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Theorem 1 Define function g(m, x) := e x TP e ' ^ ^ n ^ e ^ e ^ eas ^ integer 

such that 

I n. \ I n. \ 

> 1 - 5. 
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Observe a sequence of n interarrival times Xi, i = 1, • • • ,n. Define A = ^n n y . T/ien 

A - A 



Pr 



A 



<e\>l-5. 



Proof. Since all n inter-arrival times are i.i.d. exponential random variables with parameter A, 
the characteristic function of the average inter-arrival time, denoted by X = ^ — -, is <j>j[(t) = 



(1 - j ^x)~ n . Let Y = 2n\X = Then 4> Y = (1 - 2jt)- n which implies that F possesses a 
Chi-square distribution of degree 2n, i.e., x 2 (2n). Notice that 



A- A 



Thus 



< e 



Pr 



A x A 

< A < 



1 + e 



1-e 



2n 2n 

< F < 



1 + e 



1 - e 



A- A 



A 



< e 



2n 2n 
Pr <^ < F < 



1 + e 



1 - e 



Making use of the relation between the Chi-square distribution and the Poisson distribution (see 
[2]), we have 

Pr{F > x} = g(n,x). 
Hence, by the definition of function g (.,.), we have 



1 + e 

The proof is thus completed. 



2n 2n 
Pr < - < F < 



1 - e 



9 n , 



1 + e 



9 n 



n 



1 - e 
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Remark 1 A remarkable fact is that no a priori information about the arrival rate A is required 
to compute the sample size for pre-specified relative error bound and confidence level. 



Figure Q] shows the sample size required to obtain an estimate of the Poisson arrival rate with 
pre-specified relative accuracy (quantified by relative error bound e) and confidence level. It can 
be seen from Figure 1 that, to obtain an estimate A such that 



Pr 



A- A 



< 0.01 



0.99, 



we need to observe n = 6.6 x 10 4 inter-arrival times Xi, i = 1, ■ ■ ■ ,n and compute A = y . 
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Figure 1: Estimating Poisson Arrival Rate. Plot A corresponds to 5 = 0.05, plot B corresponds 
to 5 = 0.01, plot C corresponds to 5 = 0.001. 

4 Average Service Time 

To estimate the average of exponential service time, we have 
Theorem 2 Let e, 5 6 (0, 1). Let n be the least integer such that 

n(l + e)' 



9 \n, 



1 + 2e 



g(n,n(l + e)) >l-8. 



Let Xi, i = 1, • • • ,n be i.i.d. observations of service time. Define Jl = — -. Then 



Pr 



/' 



< e ^ > 1 - 5. 



Proof. Let Z = 2n\X = 2*2. Then Z possesses a Chi-square distribution of degree 2n, i.e., 
X 2 (2n). Notice that 



A* -A* 



/' 



< e 



Moreover, 



because 



It follows that 



V / / P 
< \i < 



1 + e 



1 + e 

2n(l + e) 
1 - e ^ 1 + 2e 

1 + e 



A 1 . . M 
< (j, < 



l-e 



< Z < 2n(l + e) 



l + 2e 



> l-e. 



Pr 



M - A* 



/' 



, 2n(l + e) 
< e !> > Pr <! / „ ' < Z < 2n 1 + e 
1 + 2e 



Making use of the relation between the Chi-square distribution and the Poisson distribution, we 
have Pr{Z > x} = g(n,x). Recall the definition of function g(., .), we have 

n(l + e)' 



pr |2n(l+ i ) <z<2n(i+e) 



1 + 2e 

The proof is thus completed 



9 n, 



l + 2e 



g(n,n(l + e)) . 
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Remark 2 It can be seen that no a priori information about the average service time \x is required 
to compute the sample size for pre- specified relative error bound and confidence level. 

Figure [2] shows the sample size required to obtain an estimate of the average service time with 
pre-specified relative accuracy (quantified by relative error bound e) and confidence level. It can 
be seen from Figure 2 that, to come up with an estimate ft such that 



Pr 



/' 



< 0.05 > > 0.95, 



we need to obtain the average of n = 1.8 x 10 3 i.i.d. observations of service time. 
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Figure 2: Estimating Average Service Time. Plot A corresponds to 5 = 0.05, plot B corresponds 
to 5 = 0.01, plot C corresponds to 5 = 0.001. 
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